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At sufficiently high densities, cold dense three-flavor quark matter is in the color-flavor locked 
(CFL) phase, in which all nine quarks pair in a particularly symmetric fashion. Once the heaviness 
of the strange quark (mass m s ) and the requirements of electric and color neutrality are taken into 
account, the CFL pattern of color superconductivity requires the pairing of quarks that would, in 
the absence of pairing, have Fermi momenta that differ by of order m 2 /n, with p the quark number 
chemical potential. This means that at sufficiently small fi, the “stress” on the pairing is large 
enough that the system can lower its energy by breaking pairs, resulting in some unconventional 
color superconductor which includes gapless excitations, spatial inhomogeneity, counter-propagating 
currents, or all three. In this paper we ask whether there is some less symmetric but still conventional 
pattern of pairing that can evade the stress. In other words, is there a pattern of pairing in which, 
once electric and color neutrality are imposed by suitable chemical potentials, pairing only occurs 
among those quarks whose Fermi momenta would be equal in the absence of pairing? We use 
graph-theoretical methods to classify 511 patterns of conventional color superconducting pairing, 
and show that none of them meet this requirement. All feel a stress, and all can be expected to 
become unstable to gapless modes at a density comparable to that at which the CFL phase becomes 
unstable. 

PACS numbers: 12.38.Mh,24.85.+p 


I. INTRODUCTION 

A. Context 

At any densities that are high enough that nucleons 
are crushed into quark matter, the quark matter that re¬ 
sults must, at sufficiently low temperatures, be in one of 
a family of color superconducting phases [I]. The essence 
of color superconductivity is quark pairing, driven by 
the BCS mechanism which operates whenever there is 
an attractive interaction between fermions at a Fermi 
surface 0- The QCD quark-quark interaction is strong 
and is attractive between quarks that are antisymmet¬ 
ric in color, so we expect cold dense quark matter to 
generically exhibit color superconductivity. Color super¬ 
conducting quark matter may well occur in the cores of 
compact stars. Except during the first few seconds af¬ 
ter their birth, compact stars are very cold compared to 
the 10’s of MeV critical temperatures expected for color 
superconductivity, and we shall therefore work at zero 
temperature throughout this paper. 

It is by now well-established that at asymptotic den¬ 
sities, where the up, down and strange quarks can be 
treated on an equal footing and the disruptive effects of 
the strange quark mass can be neglected, quark matter is 
in the color-flavor locked (CFL) phase, in which quarks 
of all three colors and all three flavors form Cooper 
pairs 111 The CFL phase is a color superconductor 
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but is an electromagnetic insulator, with zero electron 
density. To describe quark matter as may exist in the 
cores of compact stars, we need consider quark number 
chemical potentials ji of order 500 MeV at most (corre¬ 
sponding to baryon number chemical potentials fis = 3p 
of order 1500 MeV at most). The relevant range of /i 
is therefore low enough that the strange quark mass m s 
cannot be neglected: it is expected to be density depen¬ 
dent, lying between the current mass ~ 100 MeV and 
the vacuum constituent quark mass ~ 500 MeV. In bulk 
matter, as is relevant for compact stars where we are in¬ 
terested in kilometer-scale volumes, we must furthermore 
require electromagnetic and color neutrality lliil 
(possibly via mixing of oppositely charged phases) and 
allow for equilibration under the weak interactions. All 
these factors work to pull apart the Fermi momenta of the 
cross-species pairing that characterizes color-flavor lock¬ 
ing. At the highest densities, we expect CFL pairing, but 
as the density decreases the combination of nonzero m s 
and the constraints of neutrality put greater and greater 
stress on cross-species pairing, and reduce the excitation 
energies of those fermionic quasiparticles whose excita¬ 
tion would serve to ease the stress by breaking pairs H. 

If we imagine beginning in the CFL phase at asymp¬ 
totic density, reducing the density, and suppose that CFL 
pairing is disrupted by the heaviness of the strange quark 
before color superconducting quark matter is superseded 
by the hadronic phase, the CFL phase must be replaced 
by some phase of quark matter in which there is less, and 
less symmetric, pairing. The nature of this next phase 
down in density is currently not established. 

Within a spatially homogeneous ansatz, the next phase 
down in density is the gapless CFL (gCFL) phase mM 
ITl Ul ira Hi fra llfj . In this phase, quarks of all three 
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colors and all three flavors still form Cooper pairs, but 
there are regions of momentum space in which certain 
quarks do not succeed in pairing, and these regions are 
bounded by momenta at which certain fermionic quasi¬ 
particles are gapless. These gapless quasiparticles (and 
a nonzero density of electrons) make the gCFL phase an 
electric conductor. This variation on BCS pairing — in 
which the same species of fermions that pair feature gap¬ 
less quasiparticles — was first proposed for two flavor 
quark matter m and in an atomic physics context [Till . 
In all these contexts, however, the gapless paired state 
turns out in general to be unstable: it can lower its en ergy 
by the formation of counter-propagating currents 1 191 12 ( 1 . 
In the atomic physics context, the resolution of the in¬ 
stability is phase separation, into macroscopic regions of 
two phases in one of which standard BCS pairing occurs 
and in the other of which no pairing occurs 12111 . This 
picture ap pears to be supported by very recent experi¬ 
ments ( 22 , 123 . QCD context, the analogue of this 

phase separation would require coexisting phases one of 
which is positively charged and the other of which is neg¬ 
atively charged. This may in fact be the resolution to 
the instability of gapless color superconductivity in two- 
flavor quark matter |24l . In three-flavor quark matter, 
where the instability of the gCFL phase has been estab¬ 
lished in Refs. [2(i| . phase coexistence would require co¬ 
existing components with opposite color charges, in addi¬ 
tion to opposite electric charges, making it very unlikely 
that a phase separated solution can have lower energy 
than the gCFL phase mm- The gluon condensate 
phase of Ref. (23 is another possible resolution to the 
instability in two-flavor quark matter, but it requires an 
unbroken color SU( 2). In the three-flavor quark matter 
context, it seems likely that a ground state with counter- 
propagating currents is required. This could take the 
form of a crystalline color superconductor j2|| ( 23 , El] - 
the QCD analogue of a form of non-BCS pairing first con¬ 
sidered by Larkin, Ovchinnikov, Fulde and Ferrell j29j. 
Or, given that the CFL phase itself is likely augmented 
by kaon condensation (3(J, it could take the form of a 
phase in which a CFL kaon condensate carries a current 
in one direction balanced by a counter-propagating cur¬ 
rent in the opposite direction carried by gapless quark 
quasiparticles I .'I l| . This meson supercurrent phase has 
been shown to have a lower free energy than the gCFL 
phase. The investigation of crystalline color supercon¬ 
ductivity in three-flavor QCD has only just begun ]28j], 
and it remains to be seen whether such a phase can have 
a lower free energy than the meson current phase, mak¬ 
ing it a possible resolution to the gCFL instability. The 
simplest “crystal” structures do not suffice, but experi¬ 
ence in the two-flavor context [27l suggests that realistic 
crystal structures constructed from more plane waves will 
prove to be qualitatively more robust. 


B. Purpose 

The analysis of unconventional pairing in color super¬ 
conductivity is challenging, given the myriad non-BCS 
possibilities. In this context, our purpose in this paper 
is to make an exhaustive search through all the possibil¬ 
ities for conventional BCS pairing in three-flavor quark 
matter that involve pairing patterns that are less sym¬ 
metric than CFL. The first and best studied of these is 
the 2SC phase mmm, in which only quarks of two 
colors and two flavors pair. There are, however, many 
other possibilities. The question we ask in this paper is 
whether any of these patterns of pairing allow conven¬ 
tional BCS pairing to persist down to a value of the /1 
that is parametrically lower than that at which the CFL 
phase breaks down. We shall find that no such phases 
exist. 

We shall only consider patterns of pairing in which the 
Cooper pairs are antisymmetric in color, as that is the 
channel in which the interaction between quarks is attrac¬ 
tive. Furthermore, we focus on J = 0 pairing (J being the 
total angular momentum of a Cooper pair) since J = 1 
pairing between quarks of the same flavor yields gap s 
that are several orders of magnitude smaller (33 H3 |36j] . 
Antisymmetry in color, together with J = 0, requires 
antisymmetry in flavor. We shall therefore consider all 
patterns of BCS pairing that are antisymmetric in both 
color and flavor. 

The logic of our investigation is perhaps best explained 
with reference to the CFL phase itself, see Fig. ^ al- 
though the point of our investigation is to apply this 
logic to other phases. In the CFL phase, color neutral¬ 
ity requires a nonzero chemical potential /xs for the color 
charge proportional to T$ given by /ug = —to 2 /(2/x) [(J01, 
whereas electric neutrality requires /x e = 0 (37|. In the 
presence of these chemical potentials, but in the absence 
of pairing, the Fermi momenta of the nine quark species 
would be given by the effective chemical potentials plot¬ 
ted in Fig.|K The merit of this figure is that one can im¬ 
mediately see that the pattern of pairing requires pairing 
between species with differing effective chemical poten¬ 
tials, separated by 6/j, e s ~ m 2 /n. In this circumstance, 
the BCS state describing, for example, gs-bd pairing is 
constructed by first filling the gs and bd states up to a 
common Fermi surface given by the average between their 
two effective chemical potentials — which exacts a free 
energy price of order d/i 2 ff /x 2 ~ m* — and then pairing 
which gains a free energy benefit of order A 2 /x 2 . In 
the CFL phase at zero temperature the gap parameter A 
is estimated to be of order 10 to 100 MeV |3Bl- When 
the density drops sufficiently that the quark chemical po¬ 
tential /i falls below some critical value of order to 2 /A, 
then, the cost of pairing must exceed the benefit gained 
from pairing, gs-bd pairs must break, and the gCFL phase 
at that point takes over from the CFL phase. At the 
CFL-gCFL transition, however, these two phases have 
the same free energy. This means that given that the 
gCFL phase is itself unstable, there must be some other 
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FIG. 1: Schematic diagram of the effective chemical potentials 
g e a for the nine quarks as a function of the strange mass, in 
the CFL phase. We abbreviate red, green and blue by r, g 
and b and up, down and strange by u, d and s. (The effective 
chemical potentials are best thought of as the Fermi momenta 
that the quarks would have in the absence of pairing but in 
the presence of the chemical potentials g e , /X 3 and /xs that 
serve to maintain electric and color neutrality in the paired 
state.) The splitting between effective chemical potentials is 
proportional to m 2 //x. The quarks encircled with the symbol 
of the same shape (and color on line) pair with each other in 
the CFL phase. We see two instances of stressed pairing — 
i.e. pairing between quarks with differing effective chemical 
potentials. Such pairing gains a free energy benefit of order 
A 2 /x 2 but also exacts a free energy cost of order (m 2 / g) 2 g 2 = 
mi. Within a spatially homogeneous ansatz, the CFL phase 
is superseded by the gCFL phase — in which some bu-rs and 
gs-bd pairs break — when m 2 /g > 2A. However, the gCFL 
phase is itself unstable, meaning that the CFL phase must in 
fact be superseded earlier by some other phase. In this paper 
we ask whether any conventional BCS-paired state involves 
pairing only between quarks with the same effective chemical 
potential. Such a phase could take over from the CFL phase 
and persist down to a density for which m 2 /g A. We find 
no such phase. 


phase which takes over from CFL before the CFL-gCFL 
transition is reached. In this investigation, we suppose 
that this next-to-highest density phase of quark mat¬ 
ter involves conventional BCS pairing only, and involves 
pairing that is antisymmetric in color and flavor only. 
We then ask whether there is any pattern of such pairing 
in which only quarks with the same effective chemical 
potential pair. 

If the answer to this question were yes, we would then 
have a candidate pattern of pairing (or several candidate 
patterns) which could persist down to much lower densi¬ 
ties than the CFL phase, where m 2 /g^> A, without be¬ 


coming unstable towards pair-breaking, and hence with¬ 
out the need to face the question of what unconventional 
non-BCS state resolves the instability of a gapless color 
superconducting state towards the nucleation of counter- 
propagating currents. We would then have to derive and 
solve the gap equation for this stress-free pattern of pair¬ 
ing, studying it more fully than in the analysis we present 
here. 

In fact, we find no patterns of stress-free pairing. 
Among all the patterns of color and flavor antisymmetric 
BCS pairing that are possible in three-flavor quark mat¬ 
ter, stressed pairing is unavoidable. We conclude that 
even if the CFL phase were supplanted by some less sym¬ 
metric pattern of otherwise still conventional BCS pair¬ 
ing, any such phase would itself become unstable to pair 
breaking for some m 2 /g of order A. This means that un¬ 
less A /m 2 is large enough that the CFL phase persists 
down to low enough densities that it is superseded by 
hadronic matter, analysis of non-BCS options, including 
spatial inhomogeneity and/or counter-propagating cur¬ 
rents and/or gapless color superconductivity, is unavoid¬ 
able. 

The remainder of the paper is organized as follows. In 
Sec.in we present a method to implement neutrality con¬ 
ditions in superconducting matter for a generic pattern of 
pairing. This method is a generalization of that applied 
to the CFL and 2SC phases in Ref. §|. We shall assume 
that m s /g is small, but as Fig.^ m akes clear this is suffi¬ 
cient to answer the question of whether stress-free pairing 
is possible. In Sec. EH we present two simple toy models 
with two and three fermion species, respectively, in order 
to illustrate the interplay betweeen the pairing pattern 
and the neutrality conditions. In Sec. IIVI we describe the 
application of our method to three-flavor quark matter, 
including subtleties associated with imposing color neu¬ 
trality. Some of these subtleties are explained further in 
Appendix B. Sec. V contains the central argument and 
results of the paper. 


II. NEUTRALITY CONDITIONS AND 
STRESS-FREE STATES 


In this Section, we present a general form of the neu¬ 
trality conditions in a superconducting state, to leading 
order in small parameters in a sense that we shall explain. 
In subsequent sections, we shall use this formalism first 
to analyze two simple toy models and then to analyze the 
multitude of pairing patterns that are possible in quark 
matter with three colors and three flavors. 

Consider a system of N fermions with a common chem¬ 
ical potential g and K different “charge chemical poten¬ 
tials” fiQ k , meaning that the N fermions have chemical 
potentials 


K 

Vi = V + ^2<UkVQ k , i = l,...,N . (1) 

k=l 
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The coefficient qik is the Q^-charge of the *-th fermion. 
Neutrality with respect to K different Q^-charges then 
requires that the free energy f 1 satisfy 


on 

d^Q k 


= 0 , 


k = 1,..., K . 


( 2 ) 


In Sections IV and V, where we shall apply this formal¬ 
ism to quark matter, N will be 9 (quarks of three col¬ 
ors and three flavors) and K will be 3, because electric 
and color neutrality are maintained via the adjustment 
of three chemical potentials fi e , M 3 and Ms coupling to 
electric, color-T 3 , and color-Tg charges. In QCD, /j, e , M 3 
and Ms are the zeroth components of electromagnetic and 
color gauge fields, and the gauge field dynamics ensure 
that they take on values such that the matter is neutral, 
satisfying 0 {§, 23, • 

We shall only consider charges Qk which are traceless, 


'y ' Qik — 0 , (3) 

i 

as this suffices for the analysis of three-flavor quark mat¬ 
ter. If all the charges with respect to which neutrality 
must be maintained are traceless, as we assume, then 
if all the fermions have the same mass in the absence 
of pairing the matter is neutral with all p,Q k set to zero. 
We shall assume, however, that some of the fermions have 
mass m while others are massless. (In three-flavor quark 
matter m will be the strange quark mass, and we shall 
treat the light quarks as massless.) We shall always as¬ 
sume that m and all the p,Q k are much smaller than the 
common p. We shall then see that neutrality requires 
that the pQ k are parametrically of order to 2 / p. 

Our analysis can be applied to three-flavor quark mat¬ 
ter in which not all quarks pair, but it cannot be applied 
directly to two-flavor quark matter since in that context 
electric charge is not traceless. However, the question 
that we address in three-flavor quark matter in this paper 
can be answered with much less effort in the two-flavor 
case, in which Me is parametrically of order m in the ab¬ 
sence of pairing. This makes it easy to show that there 
are no stress-free patterns of pairing in two-flavor quark 
matter that are antisymmetric in color and flavor. 

We shall consider many different patterns of pairing, 
but in each case we shall simplify our notation by pre¬ 
tending that any quarks which pair do so with a common 
gap parameter A, assumed to be -C p. In an y realistic 
context, the pattern of pairing obtained by solving some 
gap equation may well have different gap parameters for 
different pairs of quarks. In such a case, the A appear¬ 
ing in our analysis should be taken to be the smallest 
nonzero gap parameter. We shall also see, momentarily, 
that the magnitude of the gap parameter A does not ac¬ 
tually play a role in our analysis; all that matters to us 
is the pattern of which quarks pair. 

The free energy of a conventional BCS state can be de¬ 
scribed as if one first assigns a common Fermi momentum 


v to all the fermions that pair with each other, and then 
reduces the free energy relative to that of this fictional 
unpaired state by a condensation energy of order A. 2 v 2 . 
If the “free energy benefit” of pairing (of order A 2 v 2 ) is 
larger than the “free energy cost” of first equalizing the 
Fermi momenta, then the BCS state under consideration 
is stable against pair-breaking. In this paper, we shall 
seek “stress-free states”, for which the pattern of pairing 
is chosen such that the free energy cost of equalizing the 
Fermi momenta is zero. 

As one can see from Fig. |TJ the CFL state has four 
sets of quarks that are “connected” by its specific pat¬ 
tern of pairing, meaning that in its analysis we will have 
to introduce four different common Fermi momenta v n . 
Anticipating the terminology from graph theory that we 
shall use in Sec. UY1 we term these sets of fermions 
components. The components of the CFL pattern are 
{fm,rs}, {bd, gs}, { rd,gu }, and {ru, gd,bs}. Unpaired 
quark matter has nine components, each consisting of a 
single quark. For the general case, we denote the compo¬ 
nents by C ni with n = 1... c where c denotes the num¬ 
ber of components. Thus, c = 4 for the CFL pattern. 
We shall assume that the ro-th component C n contains r n 
massless fermions and s n fermions with mass to. We now 
give an expression for the free energy B, which generalizes 
that in Ref. Q, assuming that fermions in components 
C n pair with gap parameter A: 


rn+Sn 


u = — X! / dpp2 r nP+S n \/p 2 + TO 2 - ^2 Mi 

77 n=l"'° V i=l 


-m 


A 2 m 2 


4-7T 2 


( 4 ) 


We have denoted the common Fermi momenta by v n \ 
these are best thought of as variational parameters and 
will be determined below by minimizing f 1 with respect 
to them. The condensation energy is written somewhat 
schematically, as it is of no consequence in our analysis as 
we now explain. As long as A <C Mj as we assume, we can 
neglect terms of order A 2 m PQ k and A 2 PQ k in the conden¬ 
sation energy, as they are much smaller than the p 2 pQ k 
terms that will turn out to determine the free energy cost 
of equalizing Fermi momenta at the values 14 . Once the 
A 2 Pq k terms are dropped, the condensation energy is as 
written, with M a number that depends on the pattern of 
pairing. (M counts the number of pairs of fermions that 
pair; Af = 12 in the CFL pattern, for example.) Upon 
making this approximation the condensation energy is 
independent of the charge chemical potentials and hence 
irrelevant in the neutrality conditions. We therefore ne¬ 
glect it henceforth. Another way of saying this is that 
our purpose is to seek patterns of pairing that minimize 
the free energy cost of equalizing the Fermi momenta of 
those quarks that pair. If we can eliminate this cost, we 
will have stress-free pairing. As we will not be evaluat¬ 
ing the free-energy benefit of pairing, the detailed form 
of the condensation energy is of no consequence. 

We now minimize U with respect to the c different com- 
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mon Fermi momenta z/„, one each for the paired fermions 
in each of the c different components C n . In so doing, we 
determine the values of the v n in terms of the chemical 
potentials and m. Because it is clear from Eq. 0 that 
the equations dfl/dv n = 0 are not coupled, we can min¬ 
imize with respect to each v n separately, yielding the 
condition 

r n +s n \ 

r n V n + Sn \/v%+m 2 - ^2 Mi I = 0 (5) 

*=I ) 


v n into the resulting expression, and differentiate with 
respect to HQ k . Both the calculation and the form of the 
result are, however, simpler if we first make the follow¬ 
ing observation. Consider U as a function of zq,..., v c 
and HQn ..., hQki where each v n itself is a function of 
/. iq 1 ,..., HQk. Then, the neutrality conditions can be 
written as 

o = AL = £“Ail + “ . (10) 

dPQ k ^ dv n dn Qk Ofj, Qk 


that must be satisfied by each of the v n . The nontrivial 
solution to this condition is 


v n 



r n +Sn 

E & 

2=1 




2 


+ {r„ 


si) m 2 . 


( 6 ) 


We now use Eq. o, and expand this expression for v n 
to linear order in the charge chemical potentials HQ k and 
to quadratic order in m. If, upon working at lowest non¬ 
trivial order in these small quantities, we were to find a 
stress-free pattern of pairing, we would then have to re¬ 
turn to this point in the analysis and investigate higher 
order corrections. As we shall show that there are no 
stress-free patterns to this order, we shall not need to 
consider higher order terms. Upon making this approxi¬ 
mation, the result for u n is most conveniently written in 
terms of “effective chemical potentials” /z® ff for each of 
the fermions which, following Ref. Q, we define as 

/i® ff = Hi for the massless fermions 

777-2 

/i® ff = Hi -for the fermions with mass m . (7) 

2 h 


We find 


r„ +s n 

E 

2=1 


( 8 ) 


We see that, to the order at which we are working, the 
common Fermi momentum v n of the fermions in the com¬ 
ponent C n is given by the arithmetic mean of the effec¬ 
tive chemical potentials of all the fermions in C n . This 
result generalizes the one of Ref. §>]. It is also straight¬ 
forward, although somewhat less instructive, to write the 
expression for u n without introducing effective chemical 
potentials: 


v n = h + 


r n +s n K 

■ E E dQ k 

2=1 k=1 


i_ rn_ 

s n 2 H 


(9) 


We now evaluate the neutrality conditions 0. The 
straightforward but lengthy way to do this is to per¬ 
form the p -integral in Eq. ©, insert the result 0 for 


The first term on the right-hand side of this equation 
vanishes as long as we are evaluating it at values of the 
v n for which dQ,/dv n = 0. Consequently, we need only 
evaluate the second term in which only the explicit HQ k ~ 
dependence of U is relevant. We see from 0 that the HQ k 
occur in U only in terms proportional to f n dp p 2 ~ z^, 
and hence conclude that the neutrality conditions take 
the form 


on 

dHQ k 



E 


V 3 n d 
3 dHQ k 


r n +s n 

E ^ ■ 

2=1 


( 11 ) 


We can simplify this upon noting from 0 that 


d 

®HQ k 


n 1*1 

E 


Hi 


rn+Sn 

'y ^ Qik 

i=1 


( 12 ) 


and by rewriting the double sum (first over components, 
then over fermions within components) as a single sum 
over all the fermions, assigning the z-th fermion a Vi given 
by the v n of the component to which it belongs. The 
resulting neutrality conditions take the form 

N 

0 = E , k=l,...,K , (13) 

i=t 


which is as one might have expected. We can further 
simplify this result by thinking of Vi as — (/x — z-'i), 
noting that our analysis is valid only to linear order in 
m 2 and HQ k an d hence, according to 0 , only to linear 
order in (/i — z/j), and therefore linearizing vf as 


v? = H 3 


3 h 2 {h - Vi) = 3 h 2 



(14) 


Upon using the traceless condition ©, the neutrality 
conditions then take on the simple form 


N 

o = E Vi qik • k = 1 ,..., K , (15) 

i=l 

In Sections IV and V we shall use the expression 0 for 
the zzj’s to apply this condition to a multitude of potential 
pairing patterns in three-flavor quark matter, after first 
applying it in a toy-model context in Section III. 
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With the notation that we have by now introduced, the 
condition for stress-free pairing can be written simply as 


Vi = nf, i = 1,... ,N . (16) 

A stress-free state has the property that the chemical po¬ 
tentials that maintain its neutrality are such that when 
we write the free energy for the BCS state we discover 
that those fermions which pair (i.e. those that lie within 
the same component) already have the same fi eS . Hence, 
the arithmetic mean in © yields simply m- Our goal, 
then, is to find a neutral, stress-free, pairing pattern, 
namely one that simultaneously fulfills Eqs. G3 and 

023- 


III. TWO TOY MODELS 


In this Section, we illustrate the interplay between the 
neutrality conditions and the pairing pattern in two toy 
models. The first model consists of two fermion species, 
one of which is massive. In the second model, we add a 
third, massless, fermion to this system in order to allow 
for more pairing patterns. 

In the first model, we have two fermion species - 
N = 2 — which means that the most general charge as¬ 
signments we can consider given the requirement of trace¬ 
lessness is simply the case with K = 1, and charge assig- 
ments qi = q and q 2 = —q. That is, we have fermions 
with chemical potentials 


= ti + q^ q , (17) 

Ii2 = fi-q/iq . (18) 


We assume the first (second) fermion to be massless 
(massive), making the effective chemical potentials 

lif = m , (19) 

= M2 - ^ • (20) 


There are only two possible patterns: pairing and no 
pairing. Without pairing, the common Fermi momenta 
are simply the effective chemical potentials, Vj = /if and 
the neutrality condition G3> yields 


flq — 



( 21 ) 


This is easy to understand: with to = 0, the system is 
neutral with fi q = 0. Switching on the mass reduces the 
effective Fermi surface of the second species reducing its 
density. Consequently, a nonzero ii q is required in order 
to restore neutrality, which here simply means equality 
of the two Fermi surfaces. The value of ii q is such that 
(if = (if = (i - m 2 /(4fi). 

When there is Cooper pairing between the two species, 
there is a single common Fermi momentum 


lif + lif 


= H ’ 


m 

If 


( 22 ) 


In this case, the neutrality condition 03, and indeed 
even the result d, are automatically satisfied, for any 
value of iiq. In this sense, this toy model is too simple, 
but let us continue with it a little farther before gener¬ 
alizing it. Since we can achieve neutrality with any ii q , 
let us choose it according to m- which makes the effec¬ 
tive potentials for the two fermions equal before pairing, 
satisfying (tTT)li and rendering the paired state stress-free. 
The fact that fi q can be varied around its stress-free value 
without upsetting neutrality is a sign that the material 
is an insulator. Varying (i q away from that in (II (ill intro¬ 
duces a stress on the pairing. There is then a free-energy 
cost to equalizing the Fermi surfaces, and the paired state 
is only favored if its condensation energy is large enough. 
The paired state will therefore be favored in a band of 
values of ii q centered on the stress-free value itTH . with 
the width of this band being controlled by the magnitude 
of A. In our analysis, we are not evaluating the conden¬ 
sation energy and so cannot evaluate the width of the 
bandgap of this insulating state. 

We can augment the toy model slightly by introducing 
“electrons”, massless fermions which by hypothesis do 
not pair and which have charge q and chemical potential 
fi q <C fi- These add a term of order to the neutral¬ 
ity condition, which is subleading relative to those we 
have already dropped in writing 03 , and so if we were 
using 03 we would not be able to evaluate the conse¬ 
quences of adding electrons. However, as we have satis¬ 
fied (1131) . in this simple model, we are justified in includ¬ 
ing electrons in the neutrality condition, and thus con¬ 
cluding that once we have introduced electrons into the 
theory we must set fi q = 0, keeping the electron density 
zero, in order to maintain neutrality. We now see that 
in this toy model augmented with the addition of “elec¬ 
trons”, the stress-free pairing condition ii q = —to 2 / fq(i) 
and neutrality condition ii q = 0 cannot both be satis¬ 
fied. Hence, pairing in neutral matter will of necessity 
be stressed, and the BCS state will only be favored if A 
is large enough or, equivalently, if for some fixed A the 
stress | iif — (if | is small enough. In fact, in this sim¬ 
ple context the quantitative criterion is easily evaluated: 
breaking a pair costs 2A in condensation energy and ben¬ 
efits | (if — (if | = m 2 /(2(i) in reduced stress, and so the 
paired state is favored as long as A > m 2 /(4/i). This toy 
model thus proves in many wavs to be analogous to the 
CFL phase, analyzed in Refs. (9, HcLIs^. 

We now consider a second toy model, in which there 
are several possible pairing patterns to consider. We add 
a third fermion species, so now N = 3, but we continue 
to keep K = 1 . We choose one of the three fermions to 
be massive, and hence consider a toy model in which the 
effective chemical potentials are 


ff 




H + qiHq 
(I + q 2 (lq 

H + q3ii g 


TO 
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components 

neutrality 

stress-free pairing 

{1},{2},{3} 

^ tf+4+i'i 

(no pairing) 

{1,2},{3} 

h q 3(73 M 

for any fj , q , if qi = q 2 

{1,3},{2} 

n — _3_i nl 

Vi 3(72 2m 

u — 1 

^ q Q3-Q1 2aa 

{2,3} ,{1} 

Ms = -3^-ir 

U — 1 mZ 

^9 ~ (73-(72 2^ 

{1,2,3} 

/iq arbitrary 

if qi = q 2 and Hi = 


TABLE I: The neutrality condition 1151 and the condition for 
stress-free pairing for a toy model with 3 fermion species. 


FIG. 2: Possible patterns of pairing in the toy model with 3 
fermion species, shown as graphs. Nonzero gap parameters Aj 
yield edges in the graphs. Below each graph, its components 
are listed. The four graphs in the second row have identi¬ 
cal components and thus are equivalent in our analysis, even 
though they would not have equivalent condensation energies. 


where 

qi + 92 + <73 = o . (24) 

Which pairing patterns are possible in this model? We 
consider the gap matrix 

/ 0 Ai A 2 \ 

M= Ai 0 A 3 , (25) 

y a 2 a 3 o ) 

where a nonzero entry in the i-th row and j-th column 
indicates that fermion i forms Cooper pairs with fermion 
j. The matrix Ai is symmetric because we are assum¬ 
ing that the Cooper pairs have zero angular momentum, 
and are thus antisymmetric in Dirac indices. We have 
written the most general possible gap matrix with diag¬ 
onal entries set to zero. In Section IV, the matrix Ai 
will be a matrix in color and flavor space, and we shall 
only consider pairing that is antisymmetric in both color 
and flavor. This means that the diagonal entries of Ai 
will vanish there, so we make this choice here too. We 
will not be solving a gap equation, and hence will only 
be interested in which of the A’s are nonzero. We shall 
therefore choose each of Ai, A 2 and A 3 to be either 
A or zero, yielding 2 3 = 8 possible patterns of pairing. 
However, some of these patterns are equivalent in their 
consequences for our analysis, as we now explain. We 
can denote each of the 8 possible patterns by a graph, 
where the vertices of the graph are the fermions and the 
edges are the nonzero gap parameters. In this simple 
example of three fermion species, this might seem nota- 
tionally superfluous but in the case of three-flavor quark 
matter, with nine fermion species, we shall find that the 
correspondence between a pairing pattern and a graph 
provides a useful tool to classify the multitude of pos¬ 
sible pairing patterns. For further explanation of the 


elementary graph-theoretical terms that we shall use in 
the following, see for instance Ref. |4Qj. All 8 graphs are 
shown in Fig. [2j and we see that the last 4 all have the 
same components — in this case, actually, the same sin¬ 
gle component. In our analysis, all that matters is the 
components since these determine which fermions have 
a common Fermi momenta. We can therefore consider 
all the graphs that have the same components as equiva¬ 
lent, even though they would have different condensation 
energies. This reduces the number of possible patterns 
that we must consider, in this example from 8 to 5. The 
number of patterns that we need consider is the num¬ 
ber of possible partitions of the set {1,2, 3,..., N} into 
disjoint, non-empty subsets. In combinatorics, this num¬ 
ber is called the Bell number Bn, and Fig. [21 shows that 
B 3 = 5. In general, Bn is given by the recursive for¬ 
mula [rTj i 

BN=J2 B *( N k 1 ) ’ (26) 

fc =0 ' ' 

along with the definition Bq = 1. This yields B± = 1, 
B 2 = 2, B 3 = 5, B 4 = 15, ... Bg = 21147, .... Return¬ 
ing to our model with N = 3, we can now discuss each 
of the 5 patterns separately, analyzing their neutrality 
conditions and looking for stress-free patterns of pairing. 
The results are collected in Table [I] 

The first row of this table describes the unpaired sce¬ 
nario. In this case, a nonzero fi q is required, as in the 
two-species toy model. 

The second, third, and fourth rows of the table describe 
patterns in which two of the fermions form Cooper pairs 
while the third remains unpaired. In the second row, for 
example, neutrality is obtained as long as 

i L q = - 2 ?-’ ( 27 ) 

\ (<?i + 92) + qi 

which yields the result quoted in the table since qi + q 2 = 
—q 3 . The analyses of the neutrality conditions in the 
third and fourth rows are similar, and in each of these 
cases a particular nonzero /i g is required. The condi¬ 
tions for stress-free pairing in the second, third and fourth 
rows of the table follow directly from m and m . For 
each of the patterns in the second, third and fourth row 
of the table, the condition that the phase is simultane¬ 
ously neutral and stress-free turns out to be that this is 



















possible only if qi = <72 = — 53 / 2 , and that it requires 
H q = —to 2 /(6qi/i). 

The fifth row of the table describes any of the patterns 
in which all three fermions pair. (This could be any of 
four different patterns, as we saw in Fig. [3) In this case, 
the system is automatically neutral, with any value of fi q . 
However, the pairing is stress-free only if qi = q 2 = — 53/2 
and fi q = —m 2 /(6qifj,). 

Our conclusion for this toy model is that if q± ^ ( 72 , 
neutral stress-free pairing is impossible, whereas if <71 = 
q 2 , we have an excess of possibilities, as all the patterns 
of pairing can be both neutral and stress free. This lat¬ 
ter result comes about because if <71 = (72 = — 93/2 and 
7 i q = — ?7Z 2 /( 6 < 7 i/u) then the unpaired phase has three 
equal effective chemical potentials. The case of three- 
flavor quark matter, to which we turn next, is more anal¬ 
ogous to the toy model with <71 7 ^ < 72 - 

IV. PAIRING PATTERNS, COLOR 
SYMMETRY, AND COLOR NEUTRALITY 

We now turn to the case of three-flavor quark matter. 
The neutrality conditions which we introduced in previ¬ 
ous sections in a general and abstract way shall now be 
those required to ensure that the quark matter is elec¬ 
trically and color neutral. We first discuss possible pair¬ 
ing patterns and the role of color and flavor symmetries. 
We then explain how color neutrality is taken into ac¬ 
count. This section provides the basis for Sec. El where 
we shall demonstrate that stress-free pairing is impossible 
in three-flavor quark matter. 


itative features of the nature of the interaction between 
quarks in QCD, we would be forced to consider a general 
9x9 gap matrix M. Restricting our attention to the 
J = 0 channel which is antisymmetric in Dirac indices 
(because J = 1 pairing yields much smaller gaps |3j j3f|) 
makes A4 a symmetric matrix specified by 45 indepen¬ 
dent gap parameters which, for our purposes, are either 
zero or nonzero. This would lead to 2 45 different gap 
matrices. However, as before many of these gap matrices 
would have the same set of components, and the number 
of patterns that have to be treated separately would be 
given by the Bell number Bg = 21147, much smaller but 
still unwieldy. 

In QCD the one-gluon exchange interaction between 
two quarks is attractive in channels which are antisym¬ 
metric under the exchange of the color of the quarks. 
We know that this is the channel in which the dominant 
attraction between quarks is to be found also by ana¬ 
lyzing the nonperturbative instanton-induced interaction 
between quarks, as well as by the qualitative argument 
that only by combining two color [3] c quarks antisymet- 
rically into a color [3] c diquark do we reduce the color 
flux. Cooper pairs that are antisymmetric in color and 
Dirac indices must also be antisymmetric in flavor in¬ 
dices, meaning that we need only consider the antisym¬ 
metric flavor channel [3]/. Therefore, we only consider 
those 9x9 gap matrices Ai in the direct product space 
[3] c 0 [3]/, which can be written as 

M = A ij Ji 0 Ij , (28) 


A. Gap matrix and color/flavor symmetries 

In quark matter with three colors and flavors, nine 
different quarks are present. Were we to ignore all qual- 


where the antisymmetric 3x3 matrices ( Ji)jk = —itijk 
and ( Ii)jk = —itijk are bases of [3] c and [3]/, respec¬ 
tively. We see that the 9x9 matrix is fully specified by 
only 9 (not 45) gap parameters, namely the 9 entries in 
Ay. Explicitly, 


0 

0 

0 

0 

-A33 

A32 

0 

A23 

— A 2 2 

0 

0 

0 

A33 

0 

—A31 

—A23 

0 

A21 

0 

0 

0 

—A32 

A31 

0 

A22 

—A21 

0 

0 

A33 

—A32 

0 

0 

0 

0 

— A13 

A12 

—A33 

0 

A31 

0 

0 

0 

A13 

0 

-An 

A32 

—A31 

0 

0 

0 

0 

— A12 

An 

0 

0 

— A23 

A22 

0 

A13 

— a 12 

0 

0 

0 

A23 

0 

— A 2 i 

— Ai 3 

0 

A n 

0 

0 

0 

—A22 

A21 

0 

A12 

—An 

0 

0 

0 

0 


with the quarks ordered ru, rd , rs, gu, gd, gs, bu, bd, 
bs. The assumptions we have made, which use only gross 
qualitative features of QCD, have excluded pairing of the 


same quarks, e.g., ru with ru , and also pairing of quarks 
of the same color, e.g., ru with rd, or with the same 
flavor, e.g., ru with bu. This is the meaning of the ze- 
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ros in the matrix A4. A pairing pattern now is uniquely 
defined by the (complex) 3x3 matrix A I;; . Therefore, 
within our method, the number of possible patterns is 
2 9 = 512, corresponding to choosing each of the nine en¬ 
tries in A ij to be either zero or A. One of these, with 
all entries in A ij set to zero, describes unpaired quark 
matter. The other 511 describe patterns of conventional 
BCS color superconducting pairing that is antisymmetric 
in both color and flavor. Note that some of these pat¬ 
terns may describe more than one phase. For example, 
both the CFL and gCFL phases would be described here 
as having a A,j with nonzero diagonal entries and zero 
off-diagonal entries, the distinction between whether the 
three diagonal entries are equal or unequal being unnec¬ 
essary to the determination of whether the pattern of 
pairing is stress-free. In our search for stress-free pairing 
described in Section V it will turn out to be most con¬ 
venient to consider all 511 patterns even though, as we 
explain now, not all are physically distinct. 

In QCD color symmetry is unbroken, but flavor sym¬ 
metry is explicitly broken both by quark masses and 
by electric charge assignments. Let us first consider, 
however, how our analysis would simplify if both flavor 
and color symmetries were unbroken. In this circum¬ 
stance, any two matrices Adi, Ad 2 that are connected by 
a rotation in color-flavor space describe the same phase. 
A color-flavor rotation U ®V where U £ SU(3) C and 
V £ SU(3)f acts on the gap matrix as 

(U®V)M = A lj UJ i U T ®VI j V T . (30) 

This rotation can equivalently be applied to the matrix 
A which transforms as 

{U®V)A = U t AV, (31) 

Hence, two matrices Ai and A 2 describe the same phase 
if there exists a rotation U®V such that Ai = U T A 2 V. 
However, for any complex matrix A 2 there is a choice of 
unitary matrices U and V such that U T A 2 V is diago¬ 
nal |42j. If flavor and color were both unbroken symme¬ 
tries we would therefore be able to restrict our attention 
to diagonal matrices A, specified by only 3 gap parame¬ 
ters An, A 2 2 and A 33 . The CFL phase corresponds to 
the case where these three gap parameters are equal . 
The gap equations for the more general case where all 
three can be nonzero and unequal have been derived and 
solved in Refs. 1GSE3 , which means that if this were 
the most general case the present analysis would not be 
necessary. 

We must, however, consider the consequences of the 
explicit breaking of flavor symmetry. In the case m u = 
m.d = 0 , m s 7 ^ 0 , which we shall consider, there appears 
to be a remaining unbroken SU (2) f symmetry. How¬ 
ever, because the u and d quarks have different electric 
charges, they will have different effective chemical po¬ 
tentials in the presence of a nonzero chemical potential 
/i e as is in general required to impose electric neutrality 
given the explicit flavor symmetry breaking introduced 


by m s 7 ^ 0. We therefore assume that there is no resid¬ 
ual SU (2) f symmetry. Then, two phases A 3 and A 2 are 
equivalent only if they are connected by a color rotation, 
i.e., if there exists a U £ SU( 3) c with Ai = U T A 2 . 
This means that we must consider both diagonal and 
non-diagonal matrices A. As we explain in the next sub¬ 
section, this complicates the analysis of the consequences 
of color neutrality. 

As an aside, note that a different analysis is required 
for the cases of superfluid 3 He and for J = 1 one- 
flavor color superconductors, where the relevant symme¬ 
try groups are SO{3)l x SO(3)s and SU(3) C x SO(3)j, 
respectively. None of these symmetries are explicitly bro¬ 
ken, and yet in both cases, non-diagonal gap matrices 
must be considered, as for example in the A phase in 
both systems |42| The reason is that although a 

general complex matrix can be diagonalized by left- and 
right- multiplication by unitary matrices, it cannot be di¬ 
agonalized if one or other matrix is required to be orthog¬ 
onal rather than unitary. In Appendix El we explicitly 
show that there is no “A phase” in the case where the 
unbroken symmetry is SU(3) C x SU(3)f. In three-flavor 
QCD, unlike in the other two cases, non-diagonal matri¬ 
ces need only be considered in the presence of explicit 
symmetry breaking. 


B. Imposing color neutrality 

In order to impose color neutrality within our ap¬ 
proach, we need to know the corresponding “color 
charge chemical potentials” HQ k , introduced in Eq. lfT]l. 
Two such chemical potentials are those originating from 
nonzero expectation values of the zeroth component of 
the gluon fields A 9 and A§. These chemical potentials 
(equivalently, gluon fields) couple to the differences be¬ 
tween color-charge densities, and so it is no surprise that 
in QCD their field equations are such as to enforce that 
they take on values corresponding to the color neutrality 
conditions ©031. However, we show in Appendix B 
that, consistent with the results of Ref. [45|, in the case of 
pairing patterns described by non-diagonal matrices A, 
color neutrality may require nonzero expectation values 
for the zeroth components of any of the gluon fields A 9 , 
a = 1,... 8 . We shall denote these by /x„, and refer to 
them as color chemical potentials, even though this ter¬ 
minology is only appropriate for /x 3 and All of them 
except H 3 and /ig appear multiplying off-diagonal Gell- 
Mann matrices in color space, meaning that they cannot 
be described as chemical potentials for suitable combina¬ 
tions of color charge densities. This in turn means that 
if any except /i 3 or /ig are nonzero, the analysis that we 
have set up in Section II is inapplicable. And, we have 
seen above that because flavor symmetry is explicitly bro¬ 
ken in QCD, we must consider non-diagonal matrices A. 
It would seem that we are in trouble. 

If we were trying to derive and solve gap equations for 
the patterns of pairing that we analyze, we would not 
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be able to evade the introduction of non-diagonal color 
“chemical potentials”. However, given our more limited 
goals, we can proceed as follows. Let A be an arbitrary 
gap matrix, in particular a non-diagonal one. For this 
gap matrix, the values of the color chemical potentials 
/i a with a = 1 ,.. .8 required for color-neutrality can be 
determined explicitly as described in Appendix H We 
set the /i a ’s to these values, and now have a color-neutral 
state. We can now apply any color rotation, after which 
we will still have a color neutral state. So, we make that 
color rotation U that diagonalizes the matrix fi a T a , where 
the T a are the Gell-Mann matrices. Since /x a T a is Her- 
mitian, such a transformation always exists. Moreover, 
since h a T a is traceless, the diagonalized matrix is also 
traceless and thus can be written as a linear combination 
of T 3 and T$, 

U/iaTaU^ =/4T 3 +/4T 8 , (32) 

where the new chemical potentials are functions of the 
old ones, H 3 = n' 3 (m ,..., /x 8 ), h’s = /4(mi, • - ■, Ms)- The 
transformation U defines a tranformation for all Gell- 
Mann matrices, T a —> UT a U\ as well as for the gap 
matrix, M —> (U ® 1) M or, equivalently, A —> U T A. 
Of course, all the new quantities depend on the chemical 
potentials Hi :..., hs- As Eq. itTTfll makes explicit, if the 
state was color neutral before this color rotation it is color 
neutral after it, assuming of course that the gap matrix 
is rotated by the same color rotation U as that applied to 
the chemical potentials. This means that for an arbitrary 
A we have found a physically equivalent order parameter 
U T A for which color neutrality requires only nonzero h 3 
and fj, 8 - 

A more formal way of explaining what we have just 
done is to note that the Cartan subalgebra of SU(3) C 
has dimension 2 , meaning that the maximal number of 
mutually commuting generators of SU( 3) c is 2 , and thus 
only two color charges can be measured simultaneously. 
This means that a phase with arbitrary pairing pattern 
A can be neutralized by the specification of only two 
linear combinations of the Ha’s. In general, these two 
combinations H 3 and /ig are non-diagonal, involving all 
the Ha s. However, by some color rotation U we find a 
new basis in which the two chemical potentials that must 
be nonzero are just h 3 and hs- Making this color rotation 
changes the pairing pattern from A to U T A. In general, 
neither A nor U T A will be diagonal. 

If we wanted to derive and solve gap equations, this re¬ 
sult would be of little practical use as in order to find the 
“simpler” state U T A explicitly one would have to have 
solved the whole problem for A. Our goals are more mod¬ 
est, however. We shall analyze the most general pairing 
pattern A, in the presence of only two color chemical po¬ 
tentials H 3 and H 81 and show that stress-free pairing is 
impossible. If the A you find interesting is one of those 
which requires other Ha’s to be nonzero, we leave it to 
you to construct U T A. Since we will have analyzed all 
possible A’s, your U T A will be equal to another A that 
we have analyzed. Thus, our analysis below suffices. 


V. STRESSED PAIRING IS UNAVOIDABLE 

In this Section, we prove that for any possible pairing 
pattern in three-flavor quark matter, there will be a mis¬ 
match in the effective chemical potentials of at least two 
fermion species that form Cooper pairs with each other. 
We end the section with a discussion of some specific pat¬ 
terns of pairing in order to illustrate our method and the 
result by example. 


A. General result 


Here we shall apply the general discussion of Section II 
to the case of interest, namely nine fermions with effective 
chemical potentials as defined in Eq. 0 given by 
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(33) 


Electric neutrality is ensured via the electric chemical 
potential /x e , while color neutrality is ensured via the 
color chemical potentials H 3 an d M 8 ; where we have nor¬ 
malized the corresponding Gell-Mann matrices as T 3 = 
diag(l/2,-1/2,0) and T 8 = diag(l/3,1/3, —2/3). Be¬ 
cause TrQ = TtT 3 = TrTg = 0, the condition 0 is 
satisfied. We consider color superconducting phases de¬ 
scribed by the gap matrix iPl). specified by the 3x3 
matrix A according to We may now proceed along 
the same lines as for the toy models in Sec. IIIII The only 
difference is the larger number of fermion species and 
the resulting larger number of possible gap matrices. Be¬ 
cause our gap matrix is specified by the nine parameters 
of A, and because we are only concerned about which pa¬ 
rameters are zero and which nonzero, we have 2 9 = 512 
patterns of pairing to consider. In order to handle this 
large number of patterns, we make use of a one-to-one 
correspondence between a gap matrix (more precisely, 
between the pattern of zero and nonzero entries in a gap 
matrix) and a graph. We can then use Maple routines 
that know the basic terms of graph theory and which, in 


11 


particular, enable us to determine the components of a 
graph. This allows us to write a program that performs 
the following steps: 

1. Create a list of all 512 pairing patterns by setting 
each entry of A either to 0 or 1. 

2. Translate each of these gap matrices into a graph 
where the vertices correspond to the nine fermions 
and the edges to nonzero entries of A. The generic 
graph (with all possible edges switched on) is shown 
in Fig. GO 

3. Determine the components of each graph. This cre¬ 
ates a list of 512 sets of componentjCi,... ,C Ci }, i = 
1,..., 512, one set for each graph and gap matrix. 
Determining the components of a graph is equiv¬ 
alent to block-diagonalizing the 9x9 matrix M 
of 1251) . For example, in the familiar CFL pattern 
in which only An, A 22 and A 33 are nonzero the 
components are Ci = {ru, gd,bs}, C 2 = {rd,gu} 7 
C 3 = {bu,rs}, and C 4 = {gs,bd}. 

4. There are many instances where several among the 
512 graphs have the same set of components. In 
fact, we find only 149 different sets of components. 
The question of whether a pattern of pairing is 
stress-free or not depends only on its set of compo¬ 
nents, so we need to analyze only 149 cases rather 
than 512. One of the sets of components contains 
only single fermions, i.e., it corresponds to unpaired 
quark matter. Thus, there are 148 different sets of 
components corresponding to the 511 patterns of 
color superconducting pairing. 

5. For each of these sets of components, compute the 
common Fermi momenta zq,..., v Ci , i = 1,..., 149, 
according to Eq. (0, and assign them to the respec¬ 
tive fermions. Here we use Eqs. CT . meaning that 
every common Fermi momentum is now a function 
of g, g e , g 3 , g s , and m s . 

6 . For each set of components, search for a simultane¬ 
ous solution (g e , g 3 , g 3 ) of the neutrality equations 
dm and the equations for stress-free pairing 113. 

The result is: For no set of components is there a triple 
(g e , g 3 , gs) that fulfills Eqs. 1151) and simultane¬ 

ously. (Except for the one that corresponds to unpaired 
quark matter. In this case, both equations are fulfilled 
trivially by ( g e , g 3 , g 8 ) = {m 2 J (4 g), 0,0).) In all 511 pos¬ 
sible patterns of BCS pairing (described by 148 different 
sets of components) if the neutrality conditions 1151) are 
satisfied there are at least two fermions which pair but 
whose values of g e s differ by of order m 2 /g. In all 511 
patterns, pairing is stressed, exacting a free energy cost 
of order (rn 2 /g) 2 g 2 = mf. Neutral stress-free BCS pair¬ 
ing is impossible. This is the central conclusion of our 
analysis. 

Before considering examples, let us review how the 
analysis of color rotations in Section IV comes into play. 



FIG. 3: Generic graph for possible pairing patterns. The 
fermions are the vertices while an edge between two vertices 
means that the corresponding fermions form Cooper pairs. A 
particular gap matrix switches off particular edges. 


Remember that for any A there is a - generally compli¬ 
cated and g a dependent - color rotated U T A for which 
color neutrality requires only nonzero g 3 and g 3 . Conse¬ 
quently, in order to discuss all possible phases it is suf¬ 
ficient to discuss all transformed matrices U T A. Any 
possible U T A, however, is included in our description, 
because all we have to know about this gap matrix is the 
corresponding set of components, and we have included 
all possible sets of components. Hence, stress-free pairing 
is impossible for any conventional color superconductor. 

In order to explain this argument somewhat differently, 
consider a specific set of components. This set may orig¬ 
inate from many different gap matrices. It will be im¬ 
possible to describe some of these gap matrices by the 
ansatz for the color chemical potentials in Eqs. ll33l) . be¬ 
cause they require nonzero g a 's other than g 3 and g 3 . 
Perhaps even all possible gap matrices corresponding to 
this specific set of components are described incorrectly. 
This would have been a problem for our method if we had 
found cases in which Eqs. m and m were simultane¬ 
ously fulfilled. Then, no direct conclusion would have 
been possible because we would not have know whether 
there was even one phase “hidden” behind the set of com¬ 
ponents which allows for stress-free pairing that is cor¬ 
rectly described using only the color chemical potentials 
g 3 and g 3 . Fortunately, our result (no stress-free pair¬ 
ing) is rigorous because the argument goes the other way 
around: any pattern of pairing is hidden behind one of 
our set of components, and since no set of components 
exhibits stress-free pairing we are able to reach a general 
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conclusion. 


B. Specific examples 

We will illustrate our result using as examples gap 
matrices A for which we know that color neutrality 
can be imposed be adjusting only p ,3 and Ms, leaving 
Mi = M 2 = M 4 = ■ • • = M 7 = 0- According to the results 
of Appendix lB 21 we may consider all matrices which have 
at most one nonzero entry per column. By choosing the 
entries of A either zero or one, this gives rise to 64 gap 
matrices (each column may either be completely zero or 
may have a 1 at three possible entries, hence there are 
4 3 = 64 options). As we prove in Appendix IQ 14 of these 
gap matrices are independent while all the others can be 
obtained from one of these 14 via a color rotation. In 
Table HI! we list these 14 matrices. 

The left half of Table EH contains all diagonal order 
parameters. One of them, here denoted by A 7 , is the 
pattern of pairing for the CFL phase. The corresponding 
graph is shown in the left panel of Fig. 0J The graphs of 
all the other diagonal order parameters in Table II are 
subgraphs of that for A 7 shown in Fig. E] i.e., they are 
obtained by removing one or more edges of this graph. 
The left half of the table includes several phases that 
are familiar from the literature. For example, the A 3 
pattern of pairing is the 2SC phase [dj|, the uSC 


and dSC phases described in Refs. [HI (HI are described 
by the patterns Ag and A 5 , and the A 2 pattern is the 
2SCus phase |H[Kj. The right half of Table ITTI shows the 
seven patterns of pairing among our 14 examples which 
have non-diagonal gap matrices. If flavor symmetry were 
unbroken, each of these patterns would be equivalent to 
one of the diagonal ones. As flavor symmetry is explicitly 
broken, however, these patterns of pairing are physically 
distinct from the ones in the left half of the table. The 
electric and color chemical potentials in the table were 
obtained from the neutrality conditions ( 113 - 
Let us discuss two of the patterns of pairing in more 
detail, explicitly showing that the results for the chemical 
potentials in Table QH do not allow for stress-free pairing. 
We do so for the patterns A 7 (CFL) and A 8 . The corre¬ 
sponding graphs are shown in Fig. 2] We take the explicit 
results from the table for /r e , /T 3 , and Ms and use them 
to evaluate the effective chemical potentials of Eqs. 63 . 
For the pattern of pairing Ag, we find (showing only the 
components with more than a single quark) 


Ci : 


eff 

M S u = M - gif 

pff 1 

Mg® =^+12% 
5m „ 

»bs = M - T^f 


tffs = M 


5 m z s 
1 ¥ 


Mb® = M ~ 6 lf 

pff 1 

Mm = A* + 12J 


(34) 

It is obvious that this state does not exhibit stress-free 
pairing. For the CFL pattern A 7 , we obtain 


J 


Ci : 


C 3 : 


n eS — n eS — n eS — 11 
M ru — Hgd. — h L bs ~ M 


6 m 


Mrf = M 


2m s 
3 u 


Me 


eff , 

ML = M + “ Me 


c 4 

r 


eff 1 hi 

Vrd = M + Me - ejf 

_ mi 

Me 6 m 


Mgu = M 


Mgs = M 
= M + 


2 m' 

“3T 


£ 

"°s 

3/i 


(35) 


Even though one linear combination of Me, M 3 and Ms is 
not fixed by the neutrality conditions, reflected here by 
the unspecified value of Me, we see that we can never¬ 
theless draw a firm conclusion. We see that the fourth 
component, C 4 , exhibits a mismatch of effective chemical 
potentials for any value of Me, meaning that stress-free 
pairing is impossible. (In a treatment that goes beyond 
the order to which we are working, the electron chemi¬ 
cal potential can be determined by including the electron 
pressure into the free energy. The result in the CFL case 
is Me = 0 urn, leading to the values for the effective 
chemical potentials illustrated in Fig.E]) There are many 
patterns of pairing among the 511 that we have analyzed 
for which, at least at the order at which we are working, 
one combination of the chemical potentials is not fixed 
by the neutrality conditions. In all such cases we have 
shown that stress-free pairing cannot arise for any value 


of the undetermined combination of chemical potentials. 


VI. CONCLUSION 

We have investigated all possible patterns of conven¬ 
tional BCS pairing in color-superconducting three-flavor 
quark matter in which the Cooper pairs are antisymmet¬ 
ric in color and flavor. We have analyzed 511 patterns of 
pairing, but it should be noted that many of them cor¬ 
respond to more than one phase, since our analysis only 
keeps track of which entries in the gap matrix are zero 
and which nonzero, not whether different nonzero entries 
are equal or unequal. Also, pushing in the opposite di¬ 
rection, not all of the 511 patterns are distinct: some are 
related by color rotation. Regardless, our analysis ap¬ 
plies to every possible pattern, only a handful of which 
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diagonal order parameters 


M3 


M8 


non-diagonal order parameters 


M3 


Ms 


Ar = 



4M 


As = 



4m 


A 2 



A 



4 M 


A3 = 



2jj 


A10 



- 


A 4 = 



2 n 


An = 



Me 


~^e + 


A 5 - 



2 [L 


2/z 


4m 


A12 = 



Me 


"Me 


2/x 


As = 



3/^ 


Me 

2 


A13 — 



2 /i e 


2/i, 


A 7 = 



Me 

2 


2/j, 


A14 = 



4m 


" 3 Ms 


Ms 


TABLE II: Order parameters for which /hi = ^2 = M 4 = ■ • • = M 7 = 0- The electric and color chemical potentials satisfying the 
neutrality constraints are given for each pattern of pairing. In those cases in which the neutrality conditions leave one linear 
combination of p e , ^3 and ps undetermined, we write the chemical potentials in terms of an unspecified fi e (or fig in the one 
case in which the undetermined combination does not involve 


have been studied before. 

At asymptotic densities, where the differences among 
the quark masses can be neglected, previous work has 
shown that the CFL phase is favored dj. As the den¬ 
sity is decreased, CFL pairing requires that the system 
pay a free energy cost that is parametrically of order 
(m^//i) 2 /i 2 = rrig, in order to gain a free energy benefit 
of order A 2 [x 2 . This guarantees that below some /i of 


order m 2 s / A, the CFL phase must break down. In this 
investigation we have shown that all 511 possible pairing 
patterns suffer the same fate. All of them exact a free 
energy cost of order ; none are stress-free. 

There are, then, three classes of possible answers to the 
question of what succeeds the CFL phase as a function 
of decreasing density. First, if m 2 /A is small enough the 
CFL phase will be succeeded directly by nuclear matter: 
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FIG. 4: Graphs for the patterns of pairing A 7 (CFL) and As, cf. Table ITTI 


the quark matter —» nuclear matter transition will oc¬ 
cur before fi drops to the critical value of order m 2 /A 
at which CFL pairing breaks down. The second pos¬ 
sibility is that the CFL phase is succeeded by a color¬ 
superconducting state with unconventional pairing, not 
simply as in BCS theory. Such a state could be one which 
is spatially inhomogeneous, one which includes quasipar¬ 
ticles which are paired but gapless, one which includes 
counterpropagating currents, one which includes meson 
condensates, or one which involves any combination of 
these possibilities. Third, it is a logical possibility that 
the CFL phase is succeeded by some among the 510 less 
symmetrically paired BCS phases that we have consid¬ 
ered, for example the 2SC phase. However, the conclu¬ 
sion of our work is that any of these must also break down 
below a /.t of order m 2 / A, meaning that any of these must 
also be succeeded by some state with non-BCS pairing. 
None of the 511 possible patterns of BCS pairing can 
survive in a regime in which /r <C m 2 / A. 
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APPENDIX A: WHY THERE IS NO “A PHASE” 
IN THE CASE OF A SU{ 3) c x SU{3) f SYMMETRY 

In this Appendix, we elaborate on the difference be¬ 
tween a three-flavor color superconductor on the one 
hand and superfluid 3 He and a one-flavor color super¬ 
conductor on the other hand. All three systems are de¬ 
scribed by an order parameter parametrized by a 3 x 3 
matrix A. The reason is that in all cases the Cooper pair¬ 
ing occurs in channels which are triplets with respect to 
two different symmetries. In a three-flavor color super¬ 
conductor these are the color and flavor (anti)triplets, in 
superfluid 3 He the angular momentum and spin triplets, 
and in a one-flavor color superconductor the color and 
total angular momentum (anti)triplets. In the latter two 
systems there is a non-diagonal order parameter matrix 

^0 0 o\ 

0 0 0, (Al) 

V 1 * 0/ 

describing the so-called A phase |Ii| |44j . This phase is 
genuinely different from all phases that can be described 
by a diagonal matrix A because there is no symmetry 
rotation that tranforms A^ into a diagonal matrix. In 
the case of 3 He, the symmetry rotations are elements of 
the group SO(3)l x SO(3)s , where SO(3)l and SO(3)s 
are the angular momentum and spin groups, respectively. 
This group acts on A as given in Eq. ED- Suppose 


v/2 
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there is a transformation U ® V with U £ SO(3)l and 

V £ SO(3)s that diagonalizes A^. This requires that 
the off-diagonal elements of the resulting matrix vanish, 

U ki Am Vtj = 0 for * 7 / j , (A2) 

u 

where Uki and vgj are the entries of the matrices U and 

V respectively. Or, explicitly, 

U 3 i(vi 2 + iv 22) = U3i(ui3 + W23) = 0 , (A3a) 

U 32 (un + iv 21) = U32(ui3 + W23) = 0 , (A3b) 

^33(^11 + iv 21) = U3 3 (ui2 + iv 22) = 0 . (A3c) 


APPENDIX B: TADPOLES AND COLOR 
NEUTRALITY 

In this Appendix we determine which color chemi¬ 
cal potentials have to be nonzero to neutralize a three- 
flavor color superconducting phase whose gap matrix A4 
is specified by OBI, with A an arbitrary 3x3 matrix. In 
QCD, these chemical potentials originate from nonzero 
expectation values of the corresponding zeroth compo¬ 
nent of the gluon fields A° a , a = 1,... 8 . In the mean-field 
approximation, the Yang-Mills equation for the gluon 
fields requires the zeroth component of the quark tad¬ 
pole contributions to vanish [33, [3f|, 


Since U is unitary and thus invertible, at least one of the 
entries U31, U32, and 1133 must be nonzero. Without loss 
of generality, we may assume that M 31 7 ! 0. Since all Vij’s 
are real, this means that v ±2 = V 22 = V 13 = V 23 = 0 . 
Therefore, 


T a ° = 0, a = 1,..., 8 . (Bl) 

Here, “quark tadpole” refers to the Feynman diagram of 
a quark loop with a gluon leg. The explicit expressions 
are given by ,39] 


VV T = 


"11 

V21V11 


V21V11 

V 21 


V31V11 

V31V21 


(A4) 


\ V 31 Vu V 31 V21 Uf-L + V%2 + V, 


33 


The orthogonality of V then implies v\ x = v\ x = 1 and 
V 21 V 11 = 0. These two conditions contradict each other. 
Consequently, there is no transformation U that di¬ 
agonalizes Aa- 

Since we only needed the v^’s to be real, the same 
proof holds for complex u t j : s, i.e., for the symmetry 
group SU(3) C x SO(3)j, relevant for a one-flavor color 
superconductor, where J is the total angular momentum. 

However, in the case of the symmetry group SU(3) C x 
SU( 3) f, hence for unitary transformations U and V, any 
matrix A can be diagonalized by U T AV, see for exam¬ 
ple Ref. |42j. In particular, a potential “A phase” order 
parameter is diagonalized by the pure flavor rotation 


U= 1 , 



/ 0 i 
0 -1 
\ V% 0 



(A5) 


where!/ £ S't/(3) c andU £ SU(3)f. With these matrices 
we find 


U T A a U 


/ 0 0 o\ 
0 0 0 

V° 0 l ) 


(A 6 ) 


Hence, if SU (3) c x SU (3)/ were unbroken in three-flavor 
QCD, the 11 A phase” would be equivalent to one with a 
diagonal A. 

As described in Section IV and manifest in Section V, 
however, the explicit breaking of SU(3)f in QCD by a 
nonzero strange quark mass means that a generic non¬ 
diagonal A is not equivalent to a diagonal A. We must 
therefore consider all 511 patterns of pairing including 
those which are analogues of the A phase. 


V J ' = -§Tr[r£W] , (B2) 

where g is the strong coupling constant. Before dis¬ 
cussing the explicit structure of the quark propagator 
S, let us point out the physical meaning of the eight 
equations m • The two tadpoles T 3 and Tg are pro¬ 
portional to color density differences, and even though 
the other tadpoles cannot be interpreted in this way this 
makes clear that satisfying the requirement that all tad¬ 
poles vanish results in a color-neutral state. One often 
replaces the precise statement that the eight g, a ’s take on 
values such that the eight 7^’s vanish by the loose state¬ 
ment that color neutrality is achieved by adjusting the 
color chemical potentials /x a such that all color densities 
vanish. As explained in Section IV, this loose statement 
can always be rendered precise by a color rotation, which 
also rotates A. 

In normal quark matter, the tadpoles vanish even with 
vanishing gluon fields A° = ... = Ag = 0. In supercon¬ 
ducting quark matter, however, nonvanishing gluon fields 
have to be introduced in order to fulfil Eqs. IBII) . Which 
gluon fields attain a nonzero expectation value depends 
on the pairing pattern. 


1. Formalism 

In Eq. IB2I) . the vertex T£ = diag(q ^T a , 
where T a are the Gell-Mann matrices, and the quark 
propagator 


5 = 




(B3) 


are given in the Nambu-Gorkov basis. We do not need 
to specify the exact form of the anomalous propagators 
iEA because they drop out in the expression for the tad¬ 
pole after taking trace over Nambu-Gorkov space. The 
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diagonal components G ± are the propagators for quarks 
and charge-conjugated quarks, 

G ± = {[G^]- 1 - ^ , (B4) 

where Gq is the tree-level quark (charge-conjugate 
quark) propagator and <f> =F the off-diagonal components 
of the quark self-energy. In general, in order to ensure 
color neutrality, all 8 color chemical potentials enter the 
(inverse) tree-level propagators, 

[G*]" 1 = 7 ± 7o (f + ha T a + He Q) - m , (B5) 

where summation over a = 1 ,..., 8 is implied, and K fJ 
is the quark four-momentum. The matrices Q and m 
are proportional to the identity in color space, and are 
diagonal in flavor space. They account for the quark elec¬ 
tric charges and masses, respectively. For the following 
qualitative discussion, we are not interested in the exact 
value of the tadpoles. We are rather interested in the 
effect of the color-flavor structure on 7),°. Thus, we may 
drop all Dirac and Nambu-Gorkov structure and write 
the tadpoles as 

T a ° = Tr jr a [F + m,T b -MHF- W T, i )" 1 M] _1 } , 

(B6) 


where F is a diagonal matrix in flavor space. We have 
replaced the (inverse) propagator by the schematic ex¬ 
pression [Gq ] -1 —» F ± n a T a and $ + — > M ., d>~ —> . 

We remark that the tadpoles remain invariant under 
color rotations of the Gell-Mann matrices and the gap 
matrix, 

T a ° [Ta , M] = r a ° [U T a u' , U M] . (B7) 

Although this is not difficult to see, we present here an 
explicit proof. The proof is basically a repeated applica¬ 
tion of the formula 

{F + UMU^)- 1 = U(F + M)~ 1 U* , (B8) 

which is valid for unitary U and U F U' = F . In the 
following, F is a matrix in flavor space, i.e., F = 1 <g> F, 
and U acts on color space, i.e., U F W = U 0 F = F; 
M can be a matrix in both color and flavor spaces, say 
M = rriij Ai ® B t with complex coefficients m, ;/ and color 
and flavor matrices A, and B 3 , respectively. Then, by 
making use of Eq. CT . 


J 


T°[UT a U\UM\ 


Tr {uT a tf [F + W UT b U * - A*jA mn UJ^ I 3 (F - UTW^ 1 UJ m U f /„] _1 } 

Tr \lJT a U t [F + m, UT b U f - A?.A mn UJ \ I 3 (F - h b T b )~ l J m U t I n ] _1 } 

Tr \UT 0 U t [F + U {^ b T b - M\F - b T b )~ l M} C/ t ]^ 1 } 

Tr |t o [F + ^ b T b -M\F-^ b T b )- l MY l ) 

T a 0 [T a ,M\, (B9) 


which is Eq. m- 


2. Evaluating tadpoles 


We may now use the general form of the tadpoles IIB(il) 
to investigate which chemical potentials are needed for 
a given gap matrix A. To this end, we have to evaluate 
the tadpoles for vanishing color chemical potentials fii = 
... = ns = 0. Let us denote these quantities by 7)^(0). 
From Eq. IIB(il) we obtain 


T a °(0) = Tr[T a (F — F _1 M) 


(BIO) 


where F = diag(/i, / 2 , /a) accounts for the explicit flavor 
symmetry breaking. A nonvanishing tadpole, '7j, 0 (0) ^ 0, 
indicates that the corresponding color chemical potential 


/r a has to be nonzero in order to obtain a color neutral 
state. The converse, namely the statement that a van¬ 
ishing tadpole, 7jj°(0) = 0 , indicates that the state can 
be color-neutralized with the corresponding color chem¬ 
ical potential /z a set to zero, is not immediate. We have 
not found a general argument that precludes the con¬ 
cern that after introducing a nonzero color chemical /x a 
that cancels the initially nonzero T)°(0), some other tad¬ 
pole 7^° could now be nonvanishing, even though 7j,°(0) 
was zero. We have, however, checked that for very many 
pairing patterns, including all those in Tableland many 
more, this concern does not arise. We have checked this 
for patterns with multiple nonzero entries per column in 
A, with complex entries in A, and for patterns with up 
to five different nonzero 7)j 0 (0)’s, which of course means 
nonzero 7jj°(0)’s that do not commute. This suggests the 
following conjecture: First, use Eq. (IbToIi to compute the 
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tadpoles with all the /j, a 's set to zero, namely the 7^(0) 
for all a = 1,... 8 . Now, keep ji a = 0 only for those 
values of a for which 7^°(0) = 0 ; for the other values 
of a, reintroduce p, a as a free variable. Then, use Eq. 
(IB 6 I) to compute all eight tadpoles as a function of 
the nonzero /j, a ’s. We conjecture that, for all pairing pat¬ 
terns A, the 7^° computed in this way will vanish for all 
values of b for which T®( 0 ) was found to vanish in the 
first step. We do not have a proof of this conjecture, but 
we know of no counterexample and, as described above, 
we have checked many examples. 

Having described the consequences of which (0) are 
zero and which nonzero, we turn now to the evaluation 
of the expression (m - We begin with the simplest pos¬ 
sible situation, setting F = 1, i.e. switching off all flavor 
symmetry breaking terms. We are left with the inversion 
of the matrix 1 — A4^A4. It can be formally inverted 
by writing the Hermitian matrix Ai'M. in terms of its 
eigenvalues Aj and the projectors onto the corresponding 
eigenspaces Vi, 


= ( Bn ) 

i=i 

where we denoted the number of different eigenvalues by 
n. Then, since {V \,..., V n } is a complete set of orthog¬ 
onal projectors, 

” 1 

(l-AfUd)- 1 = ^ —y Vi ■ (B12) 

i=l 

The projectors are given by (see for example Appendix 

A of Ref. 0) 


A - A, 


3 & 


A i A 7 


(B13) 


As a further simplification, let us assume that A4' M. has 
at most two different eigenvalues. Then, the projectors 
Vi are a linear combination of the matrices 1 and A4 
(in the case of more than two different eigenvalues, higher 
powers of M. are present). Since Tr T a = 0, we obtain 


T a °(0) oc Tr[T a M f M] = -2 Tr [T a A A 1 ] , (B14) 


where Eq. <Plt as well as Tr [T a JiJj] = — ( T a )ij and 
Tr[/,;/j] = 25ij have been used. This simple form of the 
tadpole can be applied for example to the CFL phase 
with m s = 0. In this case A = 1 and the eigenvalues of 
are 4 and 1. We immediately conclude 7^°(0) = 0 
for all a. This reflects the fact that in the absence of fla¬ 
vor symmetry breaking the CFL phase is automatically 
color neutral, i.e., there is no nonzero color chemical po¬ 
tential. 

The simplest way to introduce explicit flavor symmetry 
breaking “by hand” is to insert a flavor matrix F into Eq. 

(Em 


T a °(0) oc Tr[T a M) FM\ = - Tr[T a AF'A*] , (B15) 

where we used Tr[/ ; ; F Ij\ = Sij TrF — F tl and defined 
F' = Tr [F] ■ 1 — F. (Since the only purpose of F' is to 
break the flavor symmetry, we may drop the prime in the 
following.) The expressions for the tadpoles can now be 
evaluated explicitly: 


Ti°(0) a h Rc[A n A^] + / 2 Re[A 12 A* 2 ] + / 3 Re[A 13 A* 3 ] (B16a) 

T 2 °(0) « i (A Im[A n A 21 ] + h Im[A 12 A* 2 ] + / 3 Im[A 13 A^ 3 ]) (B16b) 

T 3 °(0) « i[/i(|A n | 2 - |A 21 | 2 ) + / 2 (|A 12 | 2 - |A 22 | 2 ) + / 3 (|A 13 | 2 - |A 23 | 2 )] (B16c) 

T 4 °(0) (x A Re[A n A 31 ] + h Re[A 12 A 32 ] + / 3 Re[A 13 A 33 ] (B16d) 

T 5 °(0) cx i (A Im[A n A 31 ] + h Im[A 12 A 32 ] + / 3 Im[A 13 A* 3 ]) (B16e) 

T 6 °(0) (x A Re[A 21 AA] + A Re[A 22 Ag 2 ] + / 3 Re[A 23 A 33 ] (B16f) 

T 7 °(0) cx i (A Im[A 21 A 31 ] + h Im[A 22 A 32 ] + f 3 Im[A 23 A* 3 ]) (B16g) 

^8°(0) ^y|[A(|Aii| 2 + |A 2 i| 2 - 2|A 3 i| 2 ) + / 2 (|Ai 2 | 2 + |A 22 | 2 - 2|A 32 | 2 ) 

+/ 3 (|A 13 | 2 + |A 23 1 2 — 2|A 33 | 2 )] . (B16h) 


Note that we have obtained Eqs. iTmTili by inserting rive them from the rigorous expression llBTTl . However 
flavor-breaking by hand, and have not been able to de- by inserting various gap matrices (including all the ex- 
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amples in Table HH and many more) explicitly into Eq. 
(IBTol) , we have checked that mm and the rigorous form 
(m yield the same conclusion. I.e., a given gap matrix 
gives rise to the same nonzero tadpoles T°( 0 ) whether 
we use fmoli or (IbTcHi . We can also use the two-flavor 
result given below in Eqs. mm as another check: Eqs. 
(TbTtI) can be obtained from Eqs. mm immediately, and 
we have been able to derive urm rigorously. 


As an example of the use of preil . we can now apply 
these equations to the CFL and gCFL phases, in both 
of which the only nonzero A’s are An, A 22 and A 33 . 
From this we see immediately that in the CFL and gCFL 
phases, the tadpoles T L °(0), X^O), T®(0),... 
vanish whereas the tadpoles T 3 °( 0 ) and ^(O) may be 
nonzero due to the explicit flavor symmetry breaking. In 
fact, in the CFL phase /1 = /2 and |An| = | A22U which 
makes T)°(0) = 0 leaving only 7?(Cl) nonzero. This is con¬ 
sistent with the result of Refs, p, 0] that the CFL phase 
is color neutral with 1 x 3 = 0 and [Xg = —to 2 /(2/z). In 
the gCFL phase, however, |An| 7 ^ |A 2 2 1 and both Tj O (0) 
and Tg’(O) are nonzero. Consequently, Eqs. IB 1GI indi¬ 
cate that in the gapless CFL phase color neutrality can 
be ensured via the introduction of suitable nonzero val¬ 
ues of 3 and /rg, as was shown explicitly in Refs. SEl 
Furthermore, Eqs. nun show tha t 1 x 1 , u?., 1 x 4 ... . 1 x 7 van¬ 
ish, as assumed in Refs. IfA lift llll ll 2 L ll. 1 l Il4 1151 Ilfl and 
demonstrated in Ref. 4: 



We can apply Eqs. unit more generally as follows. 
We note that for (a, b) being either (1,2), (4,5), or (6, 7) 
the following statement holds: both T a and % are zero if 
at least one member of each of the following three pairs 
is zero: (A m i,A n4 ), (A m 2 ,A n2 )) (A m 3 ,A n 3 ). Here, 
( m,n ) = (1,2) for (a, h) = (1,2), (to, n) = (1,3) for 
(a, h) = (4,5), and (m,n) = (2,3) for (a, h) = (6,7). 
(The argument is straightforward. Consider one of 
these pairs, say (A,„i,A„i). To make both 7^(0) and 
7)) o (0) vanish, according to Eqs. (IBl(il) we have to re¬ 
quire RefAmiA^] = ImfAmiA*!] = 0. Hence A m iA*j 
and therefore either A m i or A„i (or both) must van¬ 
ish.) Taking the case (a, b) = (1,2) as an example, we 
have learned that both Tj°( 0 ) and Tj’(O) vanish if ei¬ 
ther the first or the second (or both) entries in each col¬ 
umn of the order parameter matrix A vanish. Including 
the other two cases, we reach the following conclusion: 
Tj°(0) = T 2 °(0) = T 4 °(0) = ... = T 7 °(0) = 0 if each col¬ 
umn of the matrix A contains at most one nonzero entry. 


We can also derive an analogous result for a two-flavor 
color superconductor. This case differs from the cur¬ 
rent one in the structure of the order parameter. Since 
condensation takes place in the flavor-singlet channel, 
the order parameter is given by a three-vector A = 
(Ai,A 2 ,A 3 ) than by a 3 x 3 matrix. The gap ma¬ 
trix is Ai = A i Ji t 2 , with the antisymmetric flavor sin¬ 
glet described by the second Pauli matrix 72 . Repeat¬ 
ing the same steps as for the three-flavor case, we find 
Tr [T a M*FM] = -Tr[E] A^T a A, with F = diag(A, / 2 ), 


and thus 

7i°(0) cx (fi + / 2 ) Re[AiA 2 ] (B17a) 

T 2 °(0) (x *(/ 1 + / 2 )Im[A 1 AS] (B17b) 

^ 3 °(0) « i(A + / 2 )(|A 1 | 2 -|A 2 | 2 ) (B17 c) 

T 4 o ( 0) a (/ 1 + / 2 )Re[A 1 A*] (B17d) 

T 5 °( 0 ) a i(A + A) Im[AiA£] (B17e) 

T 6 °(0) oc (A + / 2 ) Re[A 2 Ag] (B17f) 

T 7 °(0) cx i(fi + / 2 ) Im[A 2 A 3 ] (B17g) 

T 8 °(0) oc ^±^(|A 1 | 2 + |A 2 | 2 -2|A 3 | 2 )(B17h) 

We see that T 1 °(0) = T 2 °(0) = T 4 °(0) = ... = T 7 °(0) = 0 
if the vector A contains at most one nonzero entry. 

The main results of this appendix can be summarized 
as follows: 

. The tadpoles 7?(0), 7?(0), T 4 °(0), T 5 °(0), T 6 °(0), 
7^(0) vanish if each column in the gap matrix A 
contains at most one nonzero entry. (In a two-flavor 
color superconductor, they vanish if the vector A 
contains at most one nonzero entry.) 

• If there is at least one column in the gap matrix A 
with more than one nonzero entry, the “usual” color 
chemical potentials fx 3 and fig, are not sufficient to 
ensure color neutrality. (In the two flavor case, this 
arises if the vector A has more than one nonzero 
entry. The possibility of chemical potentials other 
than fi 3 and /xg in a two-flavor color superconductor 
has been discussed in detail in Ref. |4a.1 

Note that in the two flavor case, a color rotation can al¬ 
ways be found that takes a generic gap vector A into the 
2SC pattern A = (0,0,1), meaning that chemical poten¬ 
tials other than fxg and /xg can always be avoided. In 
the three flavor case, on the other hand, a color rotation 
alone cannot turn a generic gap matrix A into one which 
has only one nonzero entry per column. We must there¬ 
fore proceed as described in Sections IV and V, making 
a color rotation that diagonalizes fx a T a but does not sim¬ 
plify A, meaning that we have to analyze an exhaustive 
list of A’s. 


APPENDIX C: SUPERCONDUCTING PHASES 
IN TABLE HU 

In this Appendix, we prove that 

1. Any of the 64 patterns of pairing described by 3 x 3 
matrices A which have at most one nonzero entry 
per column can be related to one of the 14 patterns 
of pairing in Table El by a color rotation. 

2. The 14 patterns of pairing in Table ITTI are not con¬ 
nected by any color rotation. 
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To prove the first point, we have to show that we ob¬ 
tain all 64 patterns of pairing that have at most one 
nonzero entry per column by color rotations from the 
phases in Table QH One can simply apply a color per¬ 
mutation to each of the phases. A color permutation is 
equivalent to a permutation of the rows in the matrix 
A. Therefore, Ai in Table II represents 3 matrices, and 
so do A 2 and A 3 . A 4 through A 7 represent 6 matri¬ 
ces each, etc. In total, in the same order as in Table II, 
3+3+3+6+6-1-6-1-6-1-3-1-3-1-3-1-6-1-6-1-6-1-3=63. The 64th is 
the zero matrix. 

To prove the second point, we have to show that for 
any pair (A;, Aj) chosen from among the 14 patterns of 
pairing in Table II, there is no unitary matrix U € SU (3) 
such that Aj = U Ai. This is obvious if there is one col¬ 
umn that is completely zero in A, but has a nonzero entry 
in Aj. In this case, there is no matrix — whether uni¬ 
tary or not — that can transform A; into Aj. Therefore, 
the 14 matrices can be divided into the following subsets: 
{Ai}, {A 2 }, {A 3 }, {A 4 , A§}, {A 5 , A g }, {A 6 ,Aio}, and 
{A 7 , A n , A 12 , A 13i A 14 }. For any pair (Aj, Aj) with Aj 
and Aj chosen from two different subsets, the statement 
is obviously true. The remaining cases must be treated 
separately. Suppose there is a U such that U A 4 = Ag. 
It follows that U 21 = U 22 = U 31 = U 32 = 0, where Uij are 
the entries of the matrix U. Then, 

1 |un| 1 2 3 4 5 + |wi 2 | 2 + |ui 3 | 2 U 13 U 2 3 mi 3 m 33 \ 

UU^ = U 23 UI 3 |u 23 | 2 U 2 3U*33 I ■ 

V U33U* IZ U33U23 |ti 33 | 2 J 

(Cl) 


The requirement that U be unitary leads to the contra¬ 
dictory requirements |u 23 | 2 = |u 33 | 2 = 1 and u 23 m 33 = 0. 
Consequently, A 4 and Ag are not equivalent. By symme¬ 
try, the same argument applies to the pairs (A 5 , Ag) and 
(A 6 ,Ai 0 ). Only the last of the subsets remains. Suppose 
that there is a U with U A 7 = An. This leads to a van¬ 
ishing row in U, 113 1 = it 32 = U33 = 0. Thus, U cannot 
be invertible and thus not unitary. With the same argu¬ 
ment one shows that A 7 is not equivalent to A 42 , A 43 , 
and A 44 . Now suppose that U An = A 42 . Then, there 
must be a vanishing column in U, «n = W 21 = m 3 i = 0. 
Hence this U cannot be unitary. By the same argument, 
the equivalence of An, A 43 and A 42 , A 43 is excluded. 
For U An = A 44 we find w 24 = m 22 = w 3 i = u 32 = 0. As 
shown above, this cannot be fulfilled by a unitary ma¬ 
trix. Therefore, An and A 44 are not equivalent, and, by 
symmetry, the same holds for the pairs (Ai 2 ,Ai 4 ) and 
(A 43 , A 44 ). This proves that all 14 matrices are indepen¬ 
dent. 
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